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Delta expansion and Wilson fermion in the Gross-Neveu model: 
Compatibility with linear divergence and continuum limit from inverse-mass expansion 
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We apply the 5-expansion to the Gross-Neveu model in the large N limit with Wilson fermion and investigate 
dynamical mass generation from inverse-mass expansion. The dimensionless mass M defined via the effective 
potential is employed as the expansion parameter of the bare coupling constant j3 which is partially renormalized 
by the subtraction of linear divergence. We show that 5-expansion of the 1/M series of p is compatible with 
the mass renormalization. After the confirmation of the continuum scaling of the bare coupling without fermion 
doubling, we attempt to estimate dynamical mass in the continuum limit and obtain the results converging to the 
exact value for values of Wilson parameter r G (0.8,1.0). 
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I. INTRODUCTION 

As well known, the system of real fermion fields on the lat¬ 
tice with translational invariance, chiral symmetry and locality 
has been found to contain unphysical redundancy iQl] . This no- 
go-theorem forbids any simple and straightforward fermion 
implementation into the lattice. The eldest proposal to cir¬ 
cumvent the problem is to break the chiral symmetry to the 
first order in the lattice spacing a due to the lattice inventor K. 
G. Wilson ||2l. Among other proposals iH, Wilson’s fermion 
has an advantage in the strong coupling expansion, since the 
Wilson fermion system allows an easy path to the expansion. 
The strong coupling expansion provides a simple, systematic 
and powerful computational scheme to clarify rich physics out 
of reach of continuum perturbation theory. Pertaining to the 
approach where the continuum limit is accessed by the strong 
coupling series under the crucial help of 5-expansion M, 
we like to focus on Wilson fermion system and investigate the 
recovery of the asymptotic freedom and dynamical mass gen¬ 
eration in 2D lattice Gross-Neveu model in the large N limit 

ill. 

In the formulation with auxiliary field Ux, the action of the 
lattice Gross-Neveu model reads (p = 1,2) 

S = + '4’x+^L{r + 7^)V'x 

X,fl 

-\-2ar ^ 'tpx'i/x + X '^x:'4’x'4’x 

X X 

^ X 

where i/x and g stand for the N flavor fermion on the site x 
and bare coupling constant, respectively. One choice of the 
explicit 7 matrices is 

71=0-2, 72 = 0 - 1 , 75 = 0-3 =*7172 ( 2 ) 
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where 



and tr^, = 1 (fc = 1, 2,3). The parameter r is called Wilson 
parameter and is kept non-zero to avoid fermion species dou¬ 
bling. Sm means the linearly divergent mass which is fixed 
by the one-loop computation. The fermion propagator in the 
momentum space (—Tt/a < Pfi < tt/a) reads 

Sf{p) = /I -7’ (4) 

Z-M *7/^7 sinap^ + 7 “ cosap^) 

where p takes the values 1, 2. The added extra term in 
due to Wilson behaves as (l/2)ra^^p^ nearp^ ~ 0 and is 
negligible. At a corner of Brillouin zone, p = (tt/o, tt/o) for 
example, it behaves as 4r/a and grows as a —?► 0. In addition 
of the four corners, the extra term behaves near the boundary 
of Brillouin zone as 2r/a or 4r/a and behaves as the rest mass 
which goes to infinity and decouples in the continuum limit. 

Due to the explicit breaking of the 75 symmetry at r 
0, the radiative correction for the self energy diverges lin¬ 
early and the counter term represented by Smtp'if must be ac¬ 
counted. Explicit calculation specifies that 

6m = —(2p^/a)/ (5) 

where I is given by 

j,. ^ r _(Pp __ 

(27r)2 {rX;^(l - cosp;,)}2 -hX;^sin2p^ 

( 6 ) 

By the introduction of the mass counter term, the fermion 
stays massless to all orders of perturbative expansion. 

Basic computational framework we take is the expansion 
in inverse powers of the mass suitably defined to be dimen¬ 
sionless with the combination of the lattice spacing a. Large¬ 
ness of the mass M means the largeness of the lattice spac¬ 
ing a and therefore the 1 /M expansion is equivalent with the 
strong coupling expansion. The subjects to discuss on the 
present work are as follows: First we like to show that the 







2 


mass renormalization is compatible with the large mass ex¬ 
pansion even when the 5 expansion is applied. This correct 
the wrong statement in ref. 0 that the conventional trunca¬ 
tion prescription of the ^-expansion fails to remove the linear 
divergence. We next confirm that the scaling of bare coupling 
(we do not perform full renormalization including the cou¬ 
pling constant, since our approach depends on the bare quan¬ 
tities yet) is of the physical one without fermion doubling. 
Then, we demonstrate that the mass computation can be ap¬ 
proximately carried out from large mass expansion, which is 
valid in the large lattice spacings. By such a detailed analy¬ 
sis of the (5-expansion to a solvable fermionic model, we will 
obtain one further evidence of the effectivity of ^-expansion 
combined with the large mass expansion. 


II. ESTIMATION OF THE MASS FROM 1/M EXPANSION 
A. Overview and strategy 


The mass M to be used in this work is defined through the 
effective potential V{a). In the large N limit, the fermion 
integration is of quadratic and results 


2 a^fcr — Smr 
Va^ = 




r dPp 

( 2 ^)^ 


log[ I] 


sin^ 




+ {aa + r ^ (1 - cosp^))^ 
^^ 1,2 


(7) 


smaller for larger r and larger for smaller r. For instance, at 
r = 1, M{oo, 1) = 0.46732772346 • • •. In the limit r ^ 0, 
M{oo,r) —>■ oo. Thus, 1/M expansion covers both physical 
and unphysical regions. 

Now, from (|9]l, it follows that 




1 

2 ? 


<Pp 


l-fr/MX;„(l-cosp^) 


(27r)2 {M + r X;„(l - cosp^)}^ -h Y. sin^ 


Ijr) 

M ' 


( 10 ) 


For naive fermion at r = 0, /(O) = 0 and (fTol i gives /3 = 
in \2 I ■ ii —■ Thus, d becomes a function in the 
square of M rather than M itself. In the bosonic case like 
non-linear cr models and Ising models, the inverse coupling 
constant or the inverse temperature is described in the square 
of the mass. There exists a discrepancy in the suitable mass 
parameter between the present model at r ^ 0 and bosonic 
models. 

Though unphysical region is also covered, 1 /M expansion 
for P, which we denote as /3>, is readily obtained as 


/3> 



f 1 —2r —1 + 5r^ 


( 11 ) 



The mass to be dynamically generated is given hy mjj = 
a* where a* denotes the solution of dV/da = 0. Then we 
define the dimensionless mass by 

triDa = M. ( 8 ) 

This manifests that 1 /M expansion is just an expansion effec¬ 
tive at large lattice spacings. The continuum limit is appar¬ 
ently the limit M —>• 0. 

The necessary condition of the ground state dV/da = 0 
gives the gap condition. It reads in terms of M(= aa*) as 


The above expansion becomes useless beyond the r- 
dependent convergence radius, and the small a behavior can¬ 
not be accessed. As we can see below, the (5-expansion 
changes the status in a drastic manner. 

Suppose that /3> is truncated at order n such that Pn> = 
Sfc=i The result of the (5-expansion is summarized 

by the formula 

( 12 ) 


M = aSm + 2g^ J 


dPp 


(27r)2 

M + ’'Ep(l -cosp^) 


= 25 


{M-frX;^(l-cosp^)}2-pX;^sin p^ 

dPp 


- I{r) + 


(27r)2 

M + ’'Ep(l -cosp^) 


{M-frX;„(l-cosp^)}2-pX;.sin 


(9) 


For a given positive value of g^, there corresponds one value 
of M{g^,r). For example, as 52 0 in accord with the 

asymptotic freedom, M —0 as expected. There appears, 
however, the upper limit of M in the strong coupling limit 
5 ^ ^ (X) due to the subtraction of linear divergence (In re¬ 
gion of large enough M, contribution of the counter term be¬ 
comes dominant for any non-zero r). The limit M{oo, r) is 


with the binomial coefficient 


The ^-expansion induces an order-dependent transforma¬ 
tion from Yk=ibk/M^ to the truncated series in t, 
Sfc=i bn{^)t^- Fetus use the notation D [/?„>] or simply /3„> 
for the transformation of /3. Then 

A>(() = + 

+{-l + H- \-bn ^^^f”|(14) 

Here stands for the coefficient of /3> at M“". Crucial ad¬ 
vantage of Pn> is that it exhibits the scaling behavior within 
its effective region at small t. The rigorous specification of 


n! 


k\(n — k)l 


(13) 
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FIG. 1: Plots of /? 24 > and /324< ~ l/(27r)(log(Cf) + H 2 i) at 
r = 6/5 (C = 0.4899), r = I (C = 0.5716), r = 4/5 
(C = 0.7026) and 1/2 (C = 1.2731). In the plot, we showed dashed 
line which indicates the continuum behavior at respective value of the 
Wilson parameter. We note that the harmonic constant added to the 
logarithmic term comes from the (5-expansion. 


the effective region is not known but actually, the plots shown 
in FIG. 1 exhibit the expected logarithmic continuum scal¬ 
ing of Pn>- As r decreases from 6/5, the effective region of 
Pri> grows broader. However, as found from the 4th plot at 
r = 1/2, the behavior of becomes oscillatory. The oscil¬ 


latory behavior becomes stronger for lower r. We remark that 
the oscillatory behavior is a particular property in ,5„>, the <5 
expansion of the 1/M expansion, and not observable in the 
original exact function /3(M) given by (fTOl) . 

The oscillation shows the fluctuation about the scaling be¬ 
havior and makes the estimation be too complicated. Thus, 
we understand that there are preferred values of r. Within 
the range of the preferred values, the asymptotic freedom be¬ 
havior of the bare coupling is observed and the matching of 
the behavior of Pri> with (5-expanded (3 in the scaling region 
would enable us to estimate critical quantities in the contin¬ 
uum limit. Though in ref. |@], we have wrongly stated that the 
cancellation of the first term in (fT4l i (the counter term contri¬ 
bution from the linear divergence) with the expanded series, 
the rest set of (fT4l) . is incomplete, we here correct it such that 
the renormalization of the linear divergence is effective under 
the (5-expansion . 

The most information of the scaling behavior near the con¬ 
tinuum limit is governed by the ultraviolet structure of the 
model. For example, the logarithmic behavior with the coeffi¬ 
cient l/(27r) is found from the perturbative expansion. On the 
other hand, the nonperturbative information such as the value 
of the dynamical mass cannot be reached from the sole results 
of perturbative series. In the Gross-Neveu model in the large 
N limit, the only quantity of nonperturbative nature included 
in the bare coupling is the dynamical mass to be generated. 
We like to show that the information in Pn> effective at small 
t is enough to estimate the dynamical mass rriD- To make es¬ 
timation of mu be simpler and accurate, we use a perturbative 
information of the ultraviolet divergence near the continuum 
limit in what follows. 

The behavior of the bare coupling constant is found from 
perturbative renormalization group that 

/3(M)~^log(M/C), (15) 

with unknown constant C. Since the bare coupling should 
behave as P{M) ^ log{aA[,)/{2 tt) in the M ^ 0 limit, we 
And 


niD = CAl (16) 

where A^ stands for the mass scale on the square lattice. Thus 
the estimation of the constant C directly gives the dynamical 
mass. 

The non-perturbative constant C depends on the Wilson pa¬ 
rameter and consequently too as mjy being universal. The 
non-universality of the scale Ai is natural since it depends on 
the microscopic construction of the lattice model. Analyti¬ 
cally, C is obtained by the limit, C = limM-).o exp(log M — 
2'k(3[M)). By the expansion of /3(M) given by (fTOl) in the 
mass M, we then obtain 

C = exp[27r(ci - 2c2)], (17) 


where 
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Cl 


C2 


, r r d^P 1 

lim / -—-77-T- 

M^o Ly_^ (27r)2 (M + rX;^(l - cosp^))^ + 


r d?P 


L. ( 2^)2 

- cosp^)y + y,^sIt? p^_ 


logM' 
271 - 


(18) 

(19) 


From numerical integration, we find that C = 
0.4899,0.5716,0.7026,1.2731 for r = 6/5,1,4/5,1/2, 
respectively. In experimental plots in FIG 1, we have used 
these values. 

In the large N limit, the higher order corrections to the lead¬ 
ing log is the lattice artifact. The lattice artifact contains the 
power of M and in addition the log of the mass M multiplied 
to that. These terms can be computed if we use the closed 
result (fTOl) . However, using the detailed information is not of 
real significance in our study. The aim of our study is to use 
only accessible information in the perturbation theory and the 
large M expansion. Hence, we here assume that 

P<{M) = ^\og{M/C)+R, (20) 

ZTT 

where R denotes the lattice artifact obeying limjvf-s.o 7? = 0. 
The subscript ”<” means the expansion at small M. For the 
matching of Pn> with the ^-expanded ,8<, 79[/3<] = Pn<, 
we employ the extension of the binomial coefficient by the 
Gamma functions. 


where 

/ n \ _ r(n + 1) 

V-A/ “ r(-A + l)r(n + A + l)’ 


( 21 ) 


( 22 ) 


log M as mentioned before, we simply forget it and pro¬ 
ceed in the robust manner to mimic the complicated structure 
with the simple power like corrections, 

R = cif-Pi -p C2t-P^ +■■■, (25) 


where 0 < pi < P 2 < ■ • Here we do not restrict the expo¬ 
nent pk be an integer but leave to take any positive real num¬ 
ber. In the estimation process of C via LDE, the values of 
exponents will be optimized to non-integer value for the best 
matching. 

Truncation of R to the first order gives Pn< = 
— ^{log(C't) -I- 7T„} -I- Cit~P^. The exponent of the logarith¬ 
mic term is considered as zero of double degeneracy. Thus, 
the ansatz with one-parameter obeys 


Pi 


(ilogfJ L dlogt. 


Bn< = 0. 


(26) 


In the matching region where the above LDE is valid, the 
function is also effective and approximately satisfies the 
same LDE as long as the order n is large enough. Hence, for 
large n, we deal with the same LDE for Pn>, 


0-h 


dlogt 


2 r 


Pi 


dlogt 


= 0 , 


and the integration over log t provides 


1 +Pi 


dlogf 


/3n> = -^{log(C'f) -f 77„}, 


(27) 


(28) 


Here A denotes any real number. Taking A infinitesimal in 
(EB, we obtain 1 —?> 1 and log M ^ — log t — 7T„ where the 
harmonic number 7T„ is given by 


77, 


=y-- 

Zl—/ Jc 




(23) 


The 5-expansion on /3<(M) to the order n thus provides the 
transform 


Pn<it) = -^{log(C'f) -h 77„} -h 7?„. (24) 

ZTT 

The matching of Pri> and Pn< enables us to estimate the con¬ 
stant C which directly gives the dynamical mass. The match¬ 
ing process is conveniently carried out through the use of lin¬ 
ear differential equation (LDE) to be approximately satisfied 
by Pn< ■ The construction of the LDE needs, in a strict sense, 
the information of the lattice artifact R. Though the explicit 
expansion from the gap equation (fTOl i proves the existence of 


and 


l+F/ 


dlogt. 


Pn> + ^ilogt + Hn) = -^logC. (29) 


To estimate C, we need to input values of pi and t around 
which point the LDE is considered to be satisfied. To obtain 
such optimal set of (pi,f), we employ the extension of the 
principle of minimum sensitivity (PMS) Hill. We first de¬ 
mand that the estimation should be done at the point t where 
the left-hand side of ( |29] | is stationary with respect to t. We 
further demand that the reliable estimation point should be in 
the scaling region, shown in this case as the plateau. Then, it 
is natural to employ the second derivative of the left-hand side 
of (I 29 I) be zero or approximately zero at the best estimation 
point. These conditions are written as 


1 +Pi 


d log t. 
1+Pi^ 


_l_ ^ 

277 

= 0 , 

(30) 

d 

5(2) 

Pn> 

- 0 

(31) 

dlogt. 
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The symbol in (ISTT i means the exact or the approximate 
equality (when a close point to zero exists). Note that the 
second condition OTT i is nothing but (1271) . One can first solve 
(EOll to give Pi as the function of the stationary point t, 1 /pi = 
p{t). Then, substituting the solution into the left-hand-side of 
OTb . we can obtain one or a few solutions. Among them, the 
optimal one is identified by the value of t where all relevant 
functions for (. = 0,1, 2, 3 show expected scalings. For 
example, /3„> at r = 1 shows approximate scaling about t ^ 
0.5 and optimal solution t* should be found around there. In 
these natural criteria, we can obtain only one solution at each 
order. Using the optimal solution t*, we obtain 1/p* = p{t*) 
and from ( l29l l 


C = exp 


27r(/3„> -bl/pi/3i>)|t. 


(logf* -biT„) 


(32) 


It is also possible to incorporate the next order correction 
^-P 2 xhen, LDE with which we start reads [0 -b [P 2 + 

+=0and 


d log t -I 


■(\ogt + H„) = -^log(C). 


27r 

The extended PMS conditions read 

1 r, _i d 




d 


d 


^ 271 




d 




From the first two conditions, we obtain -b p^^ = p{t) 
and (piP 2 )~^ = cr(f) as functions of t and then, from the 
third condition, optimal t = t* can be obtained within the 
observable scaling region. Following the same manner as in 
the case of one-parameter ansatz, we then obtain pj and pj 
and C from ( l33l l. 

Next order correction is difficult to incorporate, since 
the necessary higher order derivatives {£ = 6, 7) do not 
show scalings even at n = 50 which is our limit, for practical 
reason of computer facility. 


B. Estimation at r = 1 


TABLE I: Estimation result of C = 0.5716061 ■ ■ ■ in 1- and 2- 
parameter ansatz at r = 1. 


order n 20 30 40 50 

1- parameter 0.5505666 0.5610658 0.5666871 0.5688271 

2- parameter 0.5873315 0.5708970 0.5709232 
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FIG. 2: Estimation results of C = 0.5716061 ■ ■ ■ (figure (a)) and 
pi = 1 (figure (b)) in 1- and 2-parameter andatz at r = 1. Plotted 
estimation is for n = 20 to 50. 


of C-estimation as a byproduct. The limit of the sequence 
suggested is not clear yet. However, it is roughly approach¬ 
ing to the value 1, which is actually the exponent of the lead¬ 
ing order term in R: In fact, from the exact result (fTOl) . R is 
given by i? = M{const + constlogM) + 0{M^). Then 
R = const X (1/f) + 0{t~‘^), giving Pi = 1. 

We have explored the possibility of the direct estimation of 
Pi through showing scaling ^ —pi- However, we 

failed because the ratio function shows large oscillation. 


We first confine ourselves with the popular choice r = 1. 
The result of estimation up to the 2-parameter ansatz is sum¬ 
marized in Table 1 and FIG. 2(a). In FIG 2, estimation results 
of Pi^ are also plotted in (b). 

We find that the sequence of the C-estimate shows ten¬ 
dency to the exact value C = 0.5716. The speed of con¬ 
vergence is rather slow in 1-parameter ansatz. The estimate 
in 2-parameter ansatz yields accurate value but the onset of 
reliable estimation starts around order 35th. The result of pi 
estimate shown in FIG. 2(b) has been obtained from the work 


C. Estimation at r / 1 

Now, we discuss on the estimation work for r ^ 1. One 
might consider that the case r = 1/2 may provide better val¬ 
ues since the function Pn> is closer to Pn<, as seen in the last 
plot in FIG 1. However, slightly oscillates at r = 1/2 
and the derivatives would show oscillations with larger ampli¬ 
tudes. Actually, by explicit plots of (£ = 1, 2), we find 
that LDE approach does not work well due to the disturbing 



















































6 


1.04 

1.02 

1.00 

0.98 

0.96 

1.03 r 

1.02 : 
1.01 : 



I O O O 

A A 


o o oo ooo 0^0 jp-g-fj-g: 

A A X X ^ 


^0.1*4 XX' 

[i4 ^ 

; X X X X X 


20 


25 


30 


n 


40 


45 


50 


X C(1,2p)/C 

4 C(0.9,2p)/C 

o C(0.8.2p)/C 


1.00 f -*-A-i- - - j-j-,-j-*-srr 

o o o o o o 


0.99 L 


30 


35 


40 


45 


50 


1.5 

1.4 

1.3 
1.2 
1.1 
1.0 
0.9 

1.5 

1.4 
1.3 
1.2 
1.1 
1.0 
0.9 



20 


25 


30 


35 


40 


45 


50 


1 

P,(1,2p) 

Pl(0.9,2p) 

Pi(0.8,2p) 


A444J*AiSXXXxxxxx^ 

4 ^x** ^ **4Ai4i 


20 


25 


30 


n 


40 


45 


50 


FIG. 3: Estimation results of C for r = 1,0.9, 0.8 at 1- and 2- 
parameter ansatze (respectively plotted in the upper and lower). Here 
C = 0.5716 • • •, C = 0.62797 ■ • • and C = 0.7026 • • • for r = 
1, 0.9, 0.8 respectively. The plots show the ratio of estimate to the 
exact values. 


oscillation. Since the incorporation of the derivatives is cru¬ 
cial for accurate estimation, this throws us a severe problem. 
From the plots of Pn> and the derivatives, we arrive at the 
following observation. 

When r is larger than 1, the oscillation is absent but the ef¬ 
fective range of is narrow, giving less accurate estimate 
of C. When r is less than and close to 1, we observe weak 
oscillatory property in and the derivatives at low orders. 
Though this makes the estimation slightly complicated but the 
confirmation of the scaling region is possible. Around r ^ 0.8 
or slightly larger region, the scaling behavior in is roughly 

visible in low order derivatives and the estimation protocol by 
extended PMS is permitted. The value r = 0.8 is best in the 
sample three values. We confirmed that r smaller than 0.8 
makes the estimation worse due to the growing oscillation. 
This is the reason that the effective range of Wilson parameter 
is approximately found to be (0.8,1.0). The result of estima¬ 
tion is depicted in FIG 3. 

We report the result of pi estimation by pi- The result is 
shown in FIG 4. We find that, in the 1-parameter ansatz, the 
value r = 0.8 produces best estimation in the typical three 
values r = 0.8,0.9 and 1.0. At 2-parameter case, however, 
the ansatz r = 0.8 produces somewhat unstable and oscilla¬ 
tory sequence. This may be the sign that the smaller r is not 
adequate for the estimation using higher order derivatives. 

To summarize, we can say as follows: As r gets smaller, 
the region of continuum scaling observable in becomes 
wider but and its derivatives begin to show oscillation at a 


FIG. 4: Estimation results of pi = 1 for r = 1, 0.9, 0.8 at 1- and 
2-parameter ansatze (respectively plotted in the upper and lower). 
Here p\ (r, kp) denotes the estimate at Wilson parameter r with k- 
parameter ansatz. 


bit smaller value of r = 1. On the other hand, when r is larger 
(r > 1), the effective region of Pn> gets narrower and the 
scaling behavior becomes vague. We found that r G (0.8,1.0) 
provides us a good estimation up to 50th order. As in the cases 
frequently met with the Ising models, the estimation of the 
dynamical mass is better than the estimation of the exponent 
Pi- 


III. CONCLUDING REMARKS 

We first found that the cancellation of the linear divergence 
in Pn> remains effective also under the <5 expansion. As a 
consequence, the true logarithmic behavior of bare coupling 
has been observed in /3„>. We remark that the confirmation is 
explicit for the range of Wilson parameter r G (0.8,1) = I. It 
is yet unclear whether other values of r is essentially useless 
even when the order is large enough. Present 50th order study 
tells us, however, other values are not effective for practical 
use. In the range I, the estimation of the dynamical mass 
ruD is carried out in the 1- and 2-parameter ansatze and all 
the sequences indicate the convergence to the exact value. It 
is interesting to note that, from 35th to 36th orders and 43rd 
to 44th orders, rather big changes happen for r = 0.8 in 2- 
parameter ansatz. 

It would be interesting to examine the estimation with the 
use of exact value of (z = 1,2). In this case, we employ 
PMS in a looser variation of®, [l+p-i^]^W + ^^0 
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FIG. 5: Plots of estimation results for r = 1 in the 1-parameter ansatz 
with the input pi = 1 and 2-parameter ansatz with the input pi = 1 
and p 2 — 2. For comparison, the results in 1- and 2-parameter ansatz 
in the previous section with the full PMS protocol are also plotted. 


^ ~ 0 for the 2-parameter ansatz. The result for r = 1 
is shown in FIG. 5 together with the result in previous full 
PMS protocol. In the 1-parameter ansatz, the two sequences 
prove almost same accuracy (Present protocol with pi = 1 
fixed gives slightly better result). In the 2-parameter ansatz, 
the result with the exact input pi = 1 and p 2 = 2 shows bet¬ 
ter behavior to the orders of 20th or so. However, the two 
sequences tend to similar behaviors at larger orders. As this 
examination manifests itself and from the results so far ob¬ 
tained, the accuracy in the Gross-Neveu model with Wilson 
fermion is not good as in the Ising models. Actually from nu¬ 
merical tests, this is roughly understood from the behaviors of 
Pri> = a way that the scaling behavior is 

not so clear to a few tens of orders. The reason behind would 
be that non-oscillation of relevant functions needs r around 
the value r ^ 1 and in the region the lattice artifact remains 
effective. 
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